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亚 临界 雷诺 数 下 错 列 布置 三 圆柱 体 绕 流 特性 研究 


X) 5 
(广东 水 利 电力 职业 技术 学 院 ,510535 广州 ) 


摘 要 :基于 Fluent 采用 大 涡 模 拟 (LES) 方 法 ,对 亚 临界 雷诺 数 Re =3 900 条 件 下 错 列 布置 三 圆柱 
体 结 构 群 三 维 绕 流 问题 进行 了 数值 研究 。 为 验证 数值 方法 的 正确 性 ,选取 单 圆 柱 体 工 况 来 进行 相 
关 参 数 的 验证 ,发 现 结果 吻合 较 好 。 随 后 ,重点 分 析 了 间距 比 (1L/D) 此 关键 参数 对 错 列 布置 三 圆柱 
体 结 构 群 流 致 特性 的 影响 。 相 关 结 果 表 明 :在 小 间距 比 范 围 内 (1.5 达 上 六 2.0) ,中 下 游 圆柱 的 平 
均 阻 力 系数 显著 小 于 单 圆 柱 工 况 结果 ; 当 间 距 比 为 3.0 工 况 下 ,上 游 圆柱 的 时 均 压 力 系 数 分 布 与 其 
他 间距 比 工 况 结果 存在 较 大 差异 性 ;小 间距 比 工 况 下 (L/D 21.5) , 柱 体 群 尾 流 区 存在 显著 的 偏 余 
流 态 ; 三 圆柱 系统 的 边界 涡 通 量 分 布 特性 表明 , 顺 流向 涡 主 要 从 中 游 圆柱 表面 产生 ,随后 进入 流 场 
中 并 发 生变 形 与 拉 伸 。 

关键 词 :三 圆柱 ;间距 比 ; 尾 流 特 性 ;旋涡 拉 伸 与 变形 ;流体 力 系数 
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Research on flow characteristics of staggered three circular 
cylinders under subcritical Reynolds number 


LIU Qiong 


( Guangdong Polytechnic of Water Resources and Electric Engineering,510535 Guangzhou , China) 


Abstract: Based on the commercial software Fluent,the large eddy simulation ( LES) method is used to 
numerically study the three-dimensional flow around three circular cylinders in a staggered arrangement at 
Re =3 900. In order to verify the correctness of the numerical method ,a case involved in the single cylin- 
der is selected to verify the relevant parameters. Subsequently, the key parameter of the spacing ratio 
(L/D) is analyzed on the influence of the flow characteristics of the staggered three-cylinder array. Within 
the range of the small spacing ratio (1. 5 x L/D x 2.0) ,the time-averaged drag coefficient of the middle 
and downstream cylinders is smaller than that of single cylinder; when the spacing ratio is 3. 0 ,the time- 
averaged pressure coefficient distribution of the upstream cylinders 1s quite different from the results of oth- 
er cases; the formation of skewed flow in the wake of the three cylinder array can found at L/D 21.5; the 
distribution characteristics of the boundary vortex flux of the three-cylinder array are given. It shows that 


the streamwise vortex is mainly generated on the surface of the midstream cylinder and enters the flow field 
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在 实际 工程 中 三 圆柱 与 多 柱 体 结构 的 绕 流 现象 
十 分 常见 ,如 流体 流 过 海洋 管道 桥墩 结构 和 海上 平 
台 立 柱 结构 等 。 流 体 流 过 结构 产生 周期 性 流体 力 ， 
对 结构 安全 性 产生 重要 影响 。 且 多 圆柱 体 结构 群 绕 
流 之 间 的 相互 干扰 以 及 复杂 的 流动 机 制 使 得 其 比 单 
圆柱 绕 流 要 复杂 得 多 。 

目前 对 双 圆 柱 体 结构 绕 流 的 流动 特性 已 有 相关 
的 研究 成 果 , 同 样 三 圆柱 体 结构 群 绕 流 特性 的 研究 
也 有 大 量 报道 。 文 献 [1] 对 Re 2 200 时 串 列 布置 三 
同 柱 体系 统 绕 流 进行 了 数值 模拟 ,结果 表明 : 柱 体系 
统 之 间 的 间距 比 的 差异 性 可 导致 不 同 涡 旋 脱 落 模 式 
(EJERE) 。 根 据 间 距 比 ,可 将 涡 旋 脱落 模式 分 
为 了 及 类 ,并 揭示 了 串 列 布置 三 圆柱 体 之 间 的 互 扰 歼 
应 习 文 献 [2] 利 用 瞬 套 网 格 法 对 Re = 100 ,200 条 件 
下 串 列 布置 双 圆 柱 与 三 圆柱 绕 流 进行 了 数值 模拟 。 
相关 结果 表明 :圆柱 的 流体 力 系 数 和 尾 涡 发 展 模式 
受 阁 距 比 影响 较 大 ,下 游 圆柱 的 涡 劲 机 理 随 间距 比 
的 融化 与 中 上 族 圆 柱 结 采 不 同 。 

间距 比 与 雷诺 数 对 流动 形态 的 影响 、 流 体力 系 
数 和 斯 特 劳 哈 尔 数 的 影响 也 已 被 广泛 研究。 文献 
[8 对 低 雷 诺 数 下 ,间距 比 L/D =2 ® L/D =5 的 串 
列 过 圆柱 体 绕 流 的 流动 特性 进行 了 研究。 结果 表 
明 E 间 距 比 和 雷诺 数 对 多 圆柱 绕 流 的 流 态 模式 和 滴 
旋 脱 落 频 率 有 重大 影响 。 文 献 14] 使 用 开发 的 3- 
TGBS 方法 对 40 三 Re 三 160 fll L/D «4 条 件 下 并 列 三 
圆柱 结构 群 的 经 流 特 性 进行 了 研究 ,观察 到 了 8 个 
不 同 的 尾 流 模式 。 这 8 种 尾 流 模式 的 分 类 取决 于 
Re FU L/D, XL 5 ] 基于 四 步 半 隐 陈 特征 分 裂 有 限 
元 方法 (4-SICBS ) ,对 串 列 布置 三 圆柱 体 的 流 致 振动 
(FIV) 进行 了 数值 模拟 , TH] ER Eb, L/D =5,Re =80, 
120,160 时 ,在 平面 前 切 流 中 随 大 Re 的 增加 ,圆柱 在 
模 流 问 上 的 振动 啊 应 将 从 共振 现象 转变 为 涡流 干涉 。 

除了 对 三 圆柱 体 的 流动 特性 进行 研究 分 析 外 ， 
许多 学 者 还 对 四 圆柱 和 多 圆柱 阵列 的 流动 特性 也 进 
行 了 相关 人 研究。 文献 [6 对 串 列 布置 6 个 方形 圆柱 
体 群 的 绕 流 特性 进行 了 人 研究 ,分 析 间 距 比 和 雷诺 数 
等 参数 对 流 场 分 布 特 性 的 影响 ,并 阐述 了 流 态 模式 
与 间距 比 之 间 的 关系 。 文 献 [7] 研 究 了 双 圆 柱 、 三 
圆柱 和 四 圆柱 阵列 的 表面 压力 系数 变化 和 速度 功率 


谱 的 特征 。 文 献 [8 基于 二 阶 特征 分 裂 有 限 元 算法 
对 Re =100 条 件 下 六 方 柱 系统 的 统 流 特性 进行 了 二 
维 数 值 模拟 ,主要 分 析 了 间距 比 对 柱 体系 统 绕 流 特 
性 的 影响 以 及 对 流体 力 系数 的 变化 特性 的 影响 。 

现 阶段 ,对 亚 临 界 雷 诡 数 工 况 下 三 圆柱 系统 绕 
流 问 题 的 研究 仍 有 进一步 的 探索 空间 。 鉴 于 此 ,本 
人 研究 基于 商业 软件 Fluent 对 亚 临 界 雷 诺 数 (Re = 
3 900) 下 错 列 布置 三 圆柱 的 三 维 绕 流 问题 进行 了 数 
值 模拟 研究 ,主要 探究 间距 比 (LAD) 对 三 圆柱 体系 
统 的 尾 流 模 态 及 流体 力 系数 的 影响 。 此 外 ,重点 探 
讨 了 三 圆柱 的 边界 涡 通 量 分 布 特性 ,以 解释 三 圆柱 
系统 的 涡 传 输 机 制 , 涡 拉 伸 和 涡 变 形 的 机 理 。 同 时 ， 
本 人 研究 对 海洋 工程 具有 一 定 的 应 用 价值 。 在 水 流 对 
桥墩 的 冲击 作用 及 稳定 性 问题 `. 以 及 海上 钻井 平 全 
立柱 `. 运 输 管 道 等 方面 的 安全 与 稳定 性 问题 具有 重 
要 的 工程 人 赋 究 意义 。 


1 数值 计算 理论 


1.1 流体 控制 方程 


基于 大 涡 模 拟 沦 流 模型 ,通过 滤波 旺 数 将 大 尺 
度 涡 和 小 尺度 涡 分 离 , 大 尺度 涡 用 N-S 方程 直接 求 
解 , 小 尺度 涡 通 过 亚 格 子 尺 度 模 型 建立 与 大 尺度 涡 
的 关系 进行 模拟 。 经 滤波 分 析 后 ,不 可 压缩 流体 的 
Navier-Stokes 方程 为 
ð u; 


— = 0(i = 1,2,3) (1) 
ð X; 
ðu. ðu u. z Ou ðr. 
uiu MN io NE MTS UENIT ES (2) 
ot 0x; p ôx; Ox;x; Ox; 


IUP ix RR ARA b xay nsu UE SCR BE AK EIT 31 
ir x yz JT IL UE BEAT 8 S p2N Wi Ms Br SE Hs 7J sp. 为 流 
体 密度 ;iv 为 流体 黏度 ;ry = 也 — uiu; 2I AR E 
力 。 采 用 有 限 体积 法 对 控制 方程 进行 离散 ,并 采用 
SIMPLE 算 法 对 压力 与 速度 类 合 ,空间 离散 压力 项 和 
动量 项 分 别 为 二 阶 格式 和 边界 中 心 差分 格式 , 瞬 态 
方程 为 边界 二 阶 隐 式 格式 。 


1.2 无 量 纲 参 数 


相关 参数 计算 式 如 下 :阻力 系数 Co = Cn + Cn 
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704 
=2F,/ (pU HD). FH RŽ OL = Cu + Cu = 2F,/ 
(PU HD) ,压力 系数 C =2(P - P, )/ (pU) ,雷诺 数 
Re = (pUD)/v , IRE IR ARA St — fD/U, Hh, Fp 
Ej FLASH 7] 5 7173, C», Cor R28 8H 73 5 b 
性 阻力 ,Cu 与 CLOS FS283TE 73 55 8STETE7J H 和 D 分 
别 为 柱 体 轴 向 高 度 和 直径 ,P 为 静止 压力 ,P, 为 无 
穷 远 处 压力 旦 P, -0,f Aries BLUR BUS, U 为 来 
流速 度 ,p 为 流体 密度 ,> 为 流体 黏度 。 


1.3 计算 模型 与 网 格 划分 


计算 模型 中 以 x 7; p] 29 MR TC p, y. 方向 为 横流 
[I8] ,z A [6] 2 JR I9] , Je e] ea E H 2 4D,D 为 圆柱 直径 ， 
中 游 圆 柱 体 ( midstream cylinder, MC ) 圆心 坐标 设 为 
计算 域 原点 。 模 型 边界 条 件 设置 为 :入 口 边界 为 速 
ER 口 ;出 口 边界 为 压力 出 口 ;圆柱 体 结构 表面 为 无 
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2 计算 参数 选取 及 验证 


2.1 计算 参数 分 析 


以 单 圆 柱 (single cylinder, SC) THN ff] , Z5 7) 
析 了 网 格 参数 .时 间 步 长 对 斯 特务 哈 尔 数 和 平均 阻 


表 1 
Tab. 1 
算 例 时 间 步 长 /s 圆周 斑点 数 
Casel 0. 001 240 
Casel 0. 005 240 
Casel 0.010 240 
Case2 0. 005 160 
Case3 0. 005 320 
Case4 0. 005 240 
Case5 0. 005 240 
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请 移 壁面 边界 条 件 ; 其 余 周边 均 为 对 称 边 界 条 件 。 
如 图 1(a) 所 示 ,中 游 圆 柱 体 的 圆心 到 上 下 边界 的 距 
离 均 为 10D、 距 入 口 边 界 为 10D、 而 距 出 口 边界 为 
20D ,各 圆柱 体 之 间 的 间距 为 nD ,三 圆柱 的 排列 角度 
为 a=45°( 即 为 三 圆柱 中 心 点 连 线 与 外 流速 的 夹 
角 )。 计 算 域 采用 六 面体 结构 化 网 格 进行 划分 ,网 格 
尺寸 为 0.01, 展 问 市 点 数 为 41, 网 格 总 数 约 为 250 
万 ,如 图 1(b) 所 示 。 为 了 满足 y’ 的 要 求 ,通过 经 验 
公式 对 靠近 柱 体 表面 的 第 一 层 网 格 高 度 进行 计算 。 
同时 ,为 了 提高 数值 计算 的 精度 ,运用 网 格 加 密 技 术 
对 圆柱 体 周围 区 域 和 尾 流 区 域 进行 处 理 ,加 密 区 域 
为 6D x6D, 且 网 格 由 圆柱 体 癌 四 周 渐变 。 鉴 于 本 研 
究 亚 临界 雷诺 数 工 况 是 Re =3 900 ,故人 口 来 流速 度 
设 为 w=1 m/s, 流 体 密 度 p =1 000 kg/m ,流体 动力 
黏 性 系数 >=0.025 641 Pa * s, 
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(b) 网 格 划分 


计算 模型 和 网 格 划 分 


Calculation model and meshing 


力 系数 的 影响 , 详 见 表 1。 其 中 Casel 为 时 间 步 长 的 
影响 ,Case2 和 Case3 为 圆周 节点 数 的 影响 ,Case4 和 
Case5 为 边界 层 第 一 层 高 度 的 影响 。 在 计算 结果 可 
徘 的 情况 下 ,综合 考虑 数据 正确 性 与 计算 资源 的 消 
耗 性 , 本 人 研究 所 有 工 况 均 选 取 Casel 网 格 模型 与 
0. 005 s 时 间 步 长 进行 分 析 。 


Ei 
E 


时 间 步 长 .圆周 节点 和 边界 层 第 一 层 高 度 对 平均 阻力 系数 和 斯 特 劳 哈 尔 数 的 影响 


The influence of time step,circle node and first layer height of boundary layer on average drag coefficient and Strouhal number 


边界 层 第 一 层 高 度 /mm C D mean St 
0.3182 0. 994 0.21 
0.318 2 1.032 0.211 
0.318 2 1.042 0. 213 
0.3182 1. 083 0. 209 
0.3182 1:052 0. 203 
0.159 1 1.061 0.219 
1.272 8 0. 968 0. 228 
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时 间 步 长 需 符 合 库 明 数 计算 公式 要 求 :At = 
(CFL) Ax/U, 通 过 计算 可 知 时 间 步 长 取 0.005 s 可 
满足 库 明 数 要 求 。 本 人 研究 所 有 的 y' 值 均 小 于 1, 其 
y 值 可 由 公式 y =u, Ay,Zv 计算 可 得 ,其 中 尺 ,是 
壁面 摩 探 速度 ,Ay 是 第 一 层 边界 层 网 格 厚 上 度 。 为 
外 ,本 研究 数据 采集 均 在 :=290s 以 后 ,每 隔 0.005s 
采样 一 次 , 共 取 10 s。 


2.2 算 例 验证 


表 2 给 出 了 Re =3 900 时 单 融 柱 平均 阻力 系数 
和 斯 特攻 哈 尔 数 与 已 有 文献 对 比 绪 采 , RMR 
结果 与 文献 结果 吻合 较 好 ,误差 均 在 5% 以 内 ,这 论 
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证 了 本 研究 数值 方法 的 正确 性 。 图 2(a) 给 出 了 单 
圆柱 体 中 间 截 面 高 度 处 表面 的 时 均 压力 系数 与 已 有 
文献 对 比 情况 ,由 于 圆柱 表面 压力 沿 上 下 两 侧 为 对 
称 分 布 , 故 只 取 截 面 一 半 。 本 人 研究 计算 结果 与 文献 
[9-11 |] 结果 基本 吻合 ,圆周 表面 时 均 压 力 系 数 变 化 趋 
势 基本 一 致 , 仅 数值 上 存在 细小 差异 。 这 表明 本 研究 
所 采用 数值 方法 具有 一 定 的 正确 性 与 可 靠 性 。 为 外 ， 
图 2(b) 分 别 给 出 了 单 圆 柱 体 在 近 尾 流 区 (x/D = 
0. 58) 和 远 尾 流 区 (x/D =2.02) 的 时 均 速 度 分 布 情 
况 , 己 文献 [12-13 | 结 采 吻合 较 好 。 因 此 ,这 进一步 
论证 了 本 人 研究 数值 方法 的 适用 性 。 


表 2 Re=3900 工 况 下 单 圆柱 的 平均 阻力 系数 和 斯 特 劳 哈 尔 数 与 相关 文献 的 对 比 


Tab.2 Comparison of the average drag coefficient and Strouhal number of a single cylinder with related literature results at Re =3 900 


消 流 模型 Re Cp mean Cr. m. s St 
本 人 研究 结 LES 3 900 1. 032 0.262 0.211 
文献 [9|] VMS-LES 3 900 1.000 0.277 0.219 
文献 | 10 ] LES 3 900 1.016 一 0.215 
文献 [11 ] LES 3 900 1.050 E 0.210 
1:5 1.5 
We ws 
1.0 & B3 =0. 
本 研究 结果 0.0 
T o 文献 [10] -0.5 
A 文献 [1] 1.5 
1.0 EMO EAD ON CAO ONQADVON 
3 - /D—2.02 
.& 00 B 5s x 
0.0 
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0/C) 


(a) 圆柱 体 表 面 压力 系数 分 布 
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图 2 圆柱 体 表 面 压力 系数 与 尾 流 时 均 速 度 分 布 对 比分 析 


Fig.2 Comparison of the surface pressure coefficient of a cylinder and the time-averaged velocity distribution of the wake 


3 计算 结果 与 分 析 


本 章 对 排列 角度 为 45° 的 错 列 布置 三 圆柱 结构 
EEP E Yt BAE (upstream cylinder, UC ) 、 中 洲 圆 柱 
( midstream cylinder, MC) 和 和 下游 圆柱 (downstream 
cylinder, DC) 所 受 的 流体 力 及 其 周围 流 场 分 布 特性 
进行 分 析 ,并 与 单 圆 柱 工 况 结果 进行 对 比分 析 。 


3.1 流体 力 系数 分 析 
本 市 对 错 列 布置 三 圆柱 体 结构 群 的 流体 力 系数 


统计 值 随 间距 比 的 变化 规律 进行 了 分 析 , 并 与 单 圆 
柱 体 工 况 结果 进行 了 对 比 。 当 柱 体 结构 的 间距 比 小 
于 3.5 时 ,UC 的 平均 阻力 系数 呈现 逐渐 减 小 的 趋 
势 ,但 其 始终 大 于 单 圆柱 体 结果 ,如 图 3(a) 所 示 。 
同时 ,给 出 了 文献 [14] 中 错 列 双 圆 柱 的 数值 结 末 来 
对 比分 析 ; 发 现在 间距 比 L/D <2. 0 时 ,本 研究 结果 
与 文献 结果 存在 一 定 差 异 , 当 间距 比 L/D z2.5 时 ， 
文献 结果 接近 于 本 研究 结果 ,这 主要 是 由 于 文献 中 
的 计算 模型 .清流 模型 与 雷诺 数 等 参数 的 差异 所 致 。 
在 此 间距 比 范 围 内 ,MC 与 DC 的 Cnwm 值 逐渐 增 大 ， 
H L/D »2.5 时 ,MC 与 DC 的 阻力 平均 值 显著 大 
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于 单 圆 柱 体 结果 。 当 间距 比 进一步 增 大 ,三 圆柱 体 
结构 群 的 Cy,。 值 逐渐 趋 近 于 单 圆柱 体 工 况 结 

这 主要 是 由 于 圆柱 体 结构 群 之 间 的 间隙 区 增 大 ,使 
得 圆柱 体 的 尾 流 得 以 充分 发 展 ,类 似 于 单 圆柱 流动 
情形 。 值 得 注意 的 是 , 随 着 柱 体 系统 的 间距 比 变 化 ， 
UC 的 阻力 均 方 值 始终 小 于 单 圆柱 结果 ,如 图 3(b) 
所 示 。 流 体系 数 的 均 方 根 值 大 小 反映 了 相应 流体 力 
的 脉动 性 强 弱 。 由 此 可 知 ,三 圆柱 系统 中 MC 与 DC 
的 存在 对 UC 的 阻力 脉动 性 具有 抑制 作用 。 在 间距 
L/D =1.5 ~3.0 范围 内 ,MC 的 阻力 均 方 根 值 亦 小 于 
单 圆柱 结果 , 当 间 距 比 大 于 3.0 时 ,MC 的 阻力 均 方 
根 值 出 现 跳跃 性 变化 且 其 数值 显著 大 于 单 圆 柱 工 况 
结果 。 而 下 游 圆柱 的 Cj, ，. 值 在 间距 比 为 2.0 时 出 
现 跳跃 性 变化 , 且 当 间距 比 继续 增 大 时 , DC 的 
Cr IRA T. SC 结果 。UC 的 升力 系数 均 方 根 值 
随 着 间距 比 增 大 总 体 呈 现 上 升 趋势 ,在 L/D =5 T 
况 出 现 转折 点 ,如 图 3(c) 所 示 。UC 的 升力 系数 均 
刘 根 值 始终 小 于 单 圆 柱 工 况 结果 ,这 与 UC 的 阻力 
均 施 根 值 随 间距 比 变化 趋势 一 臻 。 在 间距 比 L/D = 
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12 
W 
£u. a 
- 
0.9 
© UC, 文献 [20] 


* DC, 文 献 [20] 
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(a) 平均 阻力 系数 
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(b) 阻力 均 方 根 值 


第 40 卷 
3.5 ~4.5 范围 内 ,MC 与 DC 的 升力 系数 均 方 根 值 显 
著 大 于 SC 工 况 结果 。 另 外 ,可 发 现 文献 [14 ] 的 数 
值 计算 结果 与 本 人 研究 结果 相 接 近 。 

为 进一步 分 析 柱 体系 统 的 流动 阻力 特性 ,图 4 
给 出 了 各 柱 体 的 黏 性 阻力 与 压 差 阻力 随 间 距 比 变化 
结果 ,其 中 柱 体 的 总 阻力 等 于 黏 性 阻力 与 压 差 阻力 
之 合 。 压 差 阻 力 显著 大 于 黏 性 阻力 值 旦 三 圆柱 系统 
中 柱 体 结构 的 差 压 阻力 占 主导 地 位 。UC 的 压 差 阻 
力 随 间距 比 变 化 在 L/D = 1.5 工 况 达到 最 大 值 
(C5, =1.4) , 随 间距 比 增 大 又 呈现 减少 趋势 ,如 网 
4(a) 所 示 。 同 样 ,UC 的 黏 性 阻力 在 L/D 21.5 Tht 
达到 最 小 值 , 随 间距 比 增 大 呈现 上 升 趋势 。MC 的 
压 差 阻 力 随 间距 比 变化 呈现 无 规律 形态 ,而 黏 性 阻 


力 在 间距 比 范围 ZZD=1.5 ~6.0 范围 表现 出 增 大 趋 
势 ,如 图 4(b) 所 示 。DC 的 压 差 阻 力 与 黏 性 阻力 在 
间距 比 L/D =1.5 工 况 达到 最 小 值 ( Cnm 20.78 与 
Cy, 20.027) ,如 图 4(c) 所 示 , 在 DC 阻力 成 分 分 析 
中 黏 性 阻力 数值 基本 为 正 值 , 这 与 UC .MC 的 结果 存 
在 一 定 差异 。 


13 20 25 x0 x3 10 4530 53 60 


(c) 升力 均 方 根 


图 3 错 列 布置 三 圆柱 体 流体 力 系数 随 间距 比 的 变化 


Fig.3 Variation of fluid force coefficient of three cylinders staggered arrangement with spacing ratio 
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3.2 压力 系数 分 布 


柱 体 结构 的 流体 力 系数 与 其 表面 压力 系数 存在 
相关 性 ,其 中 国 柱 表面 压力 系数 为 柱 体 的 压 差 阻力 


SC 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 6.0 
LID LD 


(b) MC (c) DC 
图 4 错 列 布置 三 圆柱 体 的 压 差 阻力 与 黏 性 阻力 随 间 距 比 的 变化 


Fig.4 Variation of pressure difference drag and viscous drag of three cylinders in staggered arrangement with spacing ratio 
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提供 主要 贡献 ,因而 分 析 三 圆柱 系统 的 柱 体 表面 压 
力 系数 具有 重要 意义 。 图 5 给 出 了 高 度 截面 Z/H = 
0.5 处 不 同 间 距 比 工 况 柱 体 表面 压力 系数 分 布 图 。 
除了 间距 比 LD=3.0 工 况 ,各 工 况 下 UC 的 压力 系 
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数 分 布 特性 相 类 似 , 如 图 5 (a) Prao TRIER TE 29 3.0 
TLF, UC 的 压力 系数 最 小 值 CC = 71.12) 7 
T SC 工 况 结果 与 其 他 间距 比 工 况 结果 ,而 其 所 对 应 
的 角度 显著 大 于 其 他 工 况 结果 ,这 说 明 。 值 得 注意 
的 是 ,小 间距 比 工 况 (LD =1.5) 的 压力 系数 最 大 值 
显著 小 于 其 他 工 况 结果 。 同 样 在 小 间距 比 工 况 , 中 
游 圆柱 的 最 大 压力 系数 显著 小 于 其 他 工 况 ,如 图 5 

1.5 1.5 


Y 
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(a) UC 


(b) 所 示 。 在 间距 比 为 3.0 时 ,MC 的 压力 系数 的 最 
小 值 显著 小 于 其 他 工 况 ,在 角度 大 于 130° 时 ,压力 系 
数 保 持 相 对 稳定 。DC 的 压力 系数 最 小 值 同样 小 于 
其 他 工 况 结果 , 且 所 对 应 的 角度 小 于 其 他 工 况 ,这 表 
明 在 此 间距 比 工 况 下 DC 的 边界 层 分 离 提 前 发 生 ， 
如 图 5(c) 所 示 。 
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图 5 不 同 工 况 下 , 错 列 布置 三 圆柱 体 的 压力 系数 分 布 特性 


Fig.5 Pressure coefficient distribution characteristics of three staggered cylinders at different cases 


3,3. 时 均 流 场 分 布 特性 

之 以 柱 体 中 间 高 度 Z/H 20.5 截面 为 分 析 对 象 ， 
术 台 将 进一步 探讨 三 种 典型 间距 比 工 况 下 (L/D = 
1GN3.0,5.0) ,三 圆柱 体 结构 群 时 均 流 场 分 布 特性 。 
在 如 间 距 比 工 况 ,三 圆柱 的 尾 流 模式 呈现 出 偏 斜 流 
模式 ,并 将 此 定义 为 偏 斜 流 态 "” ,如 图 6(a) 所 示 。 
实际 上 , 柱 体 结构 偏 斜 流 态 的 形成 与 其 自身 的 时 均 
压 考 分 布 特性 相关 ,其 中 偏 斜 流 态 的 偏 斜 度 与 偏 余 
方向 由 柱 体 结构 表面 驻 点 ( 柱 体 结构 表面 压强 最 大 
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值 与 前 应 力 为 零 所 对 应 的 位 置 ”) 位 置 决定 。 当 流 
体 流 过 柱 体 群 狭窄 的 间 际 区 后 经 过 加 速 后 形成 这 种 
偏 斜 流 态 。 小 间距 比 工 况 下 , 柱 体 之 间 的 互 扰 效应 
较为 强烈 ,其 表现 在 各 柱 体 尾 流 区 洲 涡 互相 干扰 从 
而 进一步 影响 流体 力 系 数 的 变化 特性 。UC 与 MC 
的 尾 流 区 俩 料 流 态 十 分 显著 ,而 DC BS Fé yi Dv. E 
与 长 度 显著 增 大 ,因而 DC WIA b v Ac ^E TE FÉ 
区 相对 较 远 处 ,进而 可 导致 柱 体 压 差 阻力 减少 ,这 与 
前 面 图 3(a) 的 分 析 结 来 相 类 似 。 


(a) L/D-1.5 


(b) L/D=3.0 


(c) L/D-5.0 


图 6 不 同 间距 比 工 况 下 ,三 圆柱 系统 的 时 均 涡 量 w, 和 流 线 图 分 布 


Fig.6 Time-averaged vorticity (w,) and streamline of three-cylinder system under different spacing ratios 
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当 间 距 比 增 至 3.0 时 , 柱 体 间 际 区 域 高 速 流动 
的 流体 逐渐 转变 至 正常 流速 流动 ,UC 与 MC 的 尾 流 
区 侦 料 流 态 的 俩 斜 角 度 有 所 减弱 , 且 下 洲 圆 柱 DC 
的 尾 流 区 宽度 与 长 度 缩小 至 与 MC 结果 接近 ,如 图 
6(b) 所 示 。 从 三 圆柱 系统 的 时 均 流 线 图 中 可 观察 
到 UC 的 尾 流 区 近 壁 面 附近 存在 部 分 小 尺度 的 二 次 


涡 (secondary vortex , SV) ,这 种 二 次 涡 的 存在 会 影响 


柱 体 结构 的 泄 涡 特性 。 同 样 在 MC DU Fort ECC RE Td 
处 可 显著 观察 到 小 尺度 的 二 次 涡 。 

当 间距 比 继续 增 至 5.0 时 , 柱 体 群 系统 尾 流 区 
偏 斜 流 现象 逐渐 消失 ,如 图 6(e) 所 示 。UC 的 尾 流 
形态 从 偏 斜 流 态 转变 为 正 态 , 且 其 尾 流 得 以 充分 发 
展 形成 流 涡 并 向 下 游 脱落 ,并 进一步 向 下 游 运输 。 
值得 注意 的 是 ,UC 的 尾 流 区 域 近 壁面 附近 区 域 的 小 
尺 凑 二 次 涡 现象 更 为 显著 。 同 样 MC 的 尾 流 形态 转 
变 勇 正 态 ,其 尾 涡 分 布 沿 来 流 方向 的 圆柱 中 心 线 呈 
现 对称 分 布 , 且 尾 流 区 长 度 显 著 小 于 UC 与 DC 
ac. 
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列 布置 三 圆柱 体 结构 群 的 流体 力 三 维特 性 。 图 7 给 
出 了 LAD=1.5 工 况 下 ,圆柱 体 结构 流体 力 系数 的 时 
空 分 布 。UC 的 升力 系数 与 阻力 系数 时 空 分 布 图 呈 
现 出 周期 性 变化 ,如 图 7(a) 所 示 。 这 主要 由 于 UC 
的 尾 流 充分 发 展 且 交替 脱落 ,从 而 导致 其 表面 流体 
力 呈 周期 性 变化 。 另 外 ,时 空 分 布 图 中 的 线条 颜色 
的 深浅 表示 流体 力 系数 的 数值 大 小 ,其 中 升力 系数 
时 空 分 布 图 中 细 丝 状 的 线条 宽度 与 其 自 号 的 涡 脱 落 
周期 相对 应 。 值 得 注意 的 是 ,UC 的 升 阻力 系数 时 空 
分 布 图 呈现 显 蔷 的 周期 性 变化 ,特别 是 在 沿 柱 体高 
度 方 向 ,呈现 出 剧烈 的 三 维 效 应 。MC 的 阻力 系数 
的 时 空 分 布 图 同样 呈现 出 强烈 的 三 维特 性 ,如 图 7 
(b) 所 示 , 其 升力 系数 的 时 空 分 布 表 现 出 显著 的 周 
期 性 变化 。 相 关 学 者 对 低 雷 诺 数 下 的 二 维 多 柱 
体 结 构 绕 流 特性 进行 了 数值 研究 ,但 在 柱 体 结构 的 
三 维 流动 特性 与 流体 力 特性 刻画 方面 存在 一 定局 限 
lE, DC 的 流体 力 系数 时 空 分 布 特性 呈现 出 较 大 亲 
乱 性 ,这 与 其 尾 流 区 偏 和 斜 流动 模式 密切 相关 ,如 图 
7(c) 所 示 。 
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图 7 L/D=1.5 时 ,各 柱 体 阻力 系数 (上 排 ) 和 升力 系数 (下 排 ) 时 空 分 布 特性 


Fig.7 The space-time distribution characteristics of the drag coefficient (top row) and lift coefficient (bottom row) at L/D=1.5 


3.4 ”边界 涡 通 量 分 布 特性 


由 于 圆柱 体 尾 流 特性 与 柱 体 结构 流体 力学 特性 
主要 和 尾 涡 的 演化 模式 有 关 , 因 此 ,探究 流 场 中 尾 涡 
的 来 源 十 分 有 意义 。 在 不 可 压缩 的 流体 中 , 涡 结 构 
只 能 在 无 滑 移 表 面 上 产生 并 扩散 至 流 场 中 。 边 界 高 
通 量 (boundary vorticity flux, BVF ) 表示 单位 面积 和 
时 间 在 结构 表面 上 产生 涡 旋 的 速率 ” ,由 下 式 计 算 
可 得 


o=vý n= vn Vw (3) 


AP: 》 是 结构 壁面 上 的 涡 量 梯度 张 量 ;n 是 单位 
法 回 矢 量 ;r 是 3 x1 的 矢量 。 边 界 涡 旋 通 量 的 3 个 
分 量 (o,,0,,0;) 表 示 3 WENTE (w, ,0,,0,) H 
来 源 ,其 无 量 纲 形 式 可 以 表示 为 0; 2.0 / (pU ) 

图 8 给 出 了 单 圆柱 工 况 结果 柱 体 表面 的 涡 拉 伸 
项 分 布 示 意图 。 其 中 柱 体 表 面 主 涡 ( primary vortex, 
PV ) 为 大 尺度 附着 涡 , 二 次 涡 (SV) 为 小 尺度 附着 涡 。 

图 9(a) 给 出 了 不 同 间 距 比 工 沈 下 , 错 列 布置 三 


圆柱 体 结 构 群 中 各 圆柱 体 结构 表面 的 边界 涡 通 量 分 


布 。 在 小 间距 比 工 况 下 (L/D =1.5), 涡 通 量 er 与 
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,主要 集中 于 MC 表面 ,这 表明 顺 流向 涡 与 横流 向 
涡 主要 从 MC 表面 扩散 至 流 场 中 。 边 界 涡 通 量 在 
UC 与 DC 表面 分 布 值 大 小 基本 为 0。UC 表面 涡 通 
E er. 分 布 沿 着 来 流 方 向 柱 体 中 心 线 呈现 对 称 形式 ， 
这 表明 柱 体 的 自由 剪 切 层 主要 从 柱 体 上 下 两 侧 形 成 
并 发 展 至 流 场 中 。 男 外 ,由 于 在 柱 体 群 间 隙 区 形成 
了 较 多 的 二 次 涡 , 导 人 致 涡 通 量 o. HJ) Tl 1566 273 x nib c 

当 间 距 比 增 至 3.0 时 , 涡 通 量 o, 5 o, E UC K 
面 分 布 较 弱 ,而 MC 与 DC 表面 的 涡 通 量 0 与 0, 分 
布 较 强 。 这 表明 大 尺度 的 旋涡 主要 从 中 下 游 圆 柱 表 
面 扩散 至 流 场 中 。UC 表面 涡 通 量 o, 的 分 布 较为 光 
请 平整 ,而 MC 与 DC 表面 涡 通 量 0 分布 较为 亲 乱 。 

当 涡 旋 从 圆柱 体 结构 表面 脱落 进入 流 场 时 , 涡 
旋 会 经 历 拉 伸 与 变形 阶段 ,从 而 导致 其 涡 通 量 o, 转 
换 为 o,。 因 此 ,掌握 涡 通 量 0 与 0, 之 间 的 转换 规律 
REZ, 同时 , 涡 通 量 o, ATI JE IM h A EE 


用 3 为 了 进一步 分 析 涡 运输 机 理 ,引入 涡 拉 伸 项 
PE. o, CoU/0:0 和 涡 变 形 项 P2 = o, aUa) + 


oUa 用 于 研究 顺 流向 涡 co, fria LSU LoT 
算 会 式 为 

| 2 _ _ 

Ardo., = _ 

C dt ox — ^ y 0z 


ON 图 9(b) 给 出 了 错 列 布置 三 圆柱 体 结构 群 表 面 
PE 和 P2 分 布 特性 。 随 着 间距 比 变 化 ,三 圆柱 体 结 
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(a) 边界 涡 通 量 分 布 
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构 群 涡 拉 伸 和 变形 机 理 具有 较 大 的 差异 。 实 际 上 ， 
P1 和 P2 项 分 别 反 映 顺 流 癌 涡 w, 的 拉 伸 程度 与 变 
形 程度 ,Pl + P2 表示 w, 的 拉 伸 项 与 变形 项 综合 影 
啊 程 度 。P1 与 P2 项 的 分 布 沿 着 三 圆柱 表面 呈现 反 
对 称 分 布 。 在 小 间距 比 工 况 下 ,Pl 与 P2 项 主 涡 
(PV) 分布 主要 集中 于 MC 与 DC 表面 。P1 + P2 合 
项 的 改变 会 导致 主 涡 分 布 强度 减弱 ,其 中 主要 表现 
为 涡 拉 伸 机 制 与 涡 变 形 机 制 相 互 抵 消 , 最 终 诱 导 SV 
的 形成 。 而 Pl1 + P2 合 项 的 PV 分 布 在 圆柱 体 表面 
几乎 为 0。 当 间距 比 增 至 3.0 时 ,三 圆柱 表面 的 P1 
与 P2 项 PV 分 布 较 多 ;不 同 的 是 , MC 与 DC 表面 的 
大 尺度 P1 与 P2 项 PV 与 SV 分 布 均 较 多 。 三 圆柱 
表面 的 Pl + P2 合 项 主 涡 分 布 与 二 次 涡 分 布 基本 
为 0。 


PV 


图 8 单 圆柱 工 况 涡 拉 伸 项 分 布 示 意图 


Fig.8 The distribution of the vortex stretching 


term of the circular cylinder 
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(b) 涡 拉 伸 与 涡 变 形 分 布 
图 9 当 ZD=15 与 3.0 时 ,三 圆柱 系统 的 边界 润 流通 量 瞬 时 分 布 以 及 润 拉 伸 项 和 涡 变 形 项 分 布 


Fig.9 The instantaneous distribution of the boundary vortex flux and the distribution of the vortex stretching and vortex turning 


terms of the three — cylinder array at L/D 21.5 and 3.0 


4 £j ie 


本 研究 对 Re =3 900 下 ,排列 角度 为 a 2 45? fft 
置 三 圆柱 体 结构 群 的 三 维 绕 流 问题 进行 了 数值 计 
算 。 分 析 了 三 圆柱 系统 的 流体 力 系数 随 间距 比 的 分 


布 特性 ,随后 探究 了 流 场 分 布 特性 ,并 解释 了 相关 涡 
动机 理 ,重点 揭示 了 涡 旋 发 展 传播 的 演化 规律 。 主 
要 结论 如 下 。 

1) 中 下 游 加 柱 体 的 平均 阻力 系数 随 间距 比 变 化 
先 呈 现 上 升 趋势 ,随后 逐渐 接近 单 圆柱 体 工 况 结 
而 上 游 贺 柱 的 阻力 系数 均 方 根 值 与 升力 系数 均 方 根 
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值 随 间距 比 变化 始终 小 于 单 加 柱 体 工 况 结果 。 间 距 
比 为 3.0 工 沉 的 时 均 压 力 系 数 分 布 模式 与 其 他 工 况 
结果 存在 较 大 差 寞 性 。 

2) 在 小 间距 比 工 况 ,三 圆柱 系统 的 尾 流 区 存在 
显 若 的 俩 料 流 态 模式 , 随 间距 比 增 大 ,这 种 俩 斜 流 态 
模式 逐渐 转变 正 第 态 。 同 样 在 小 间距 比 工 况 下 , 柱 
体 间 际 区 存在 显著 的 小 尺度 二 次 涡 , 对 柱 体 之 间 的 
互 扰 效应 具有 重要 影响 。 

3) 三 圆柱 表面 的 边界 消 通 量 分 布 模式 与 其 旋涡 
来 源 以 及 运输 模式 存在 紧密 联系 。 在 小 间距 比 工 况 
下 , 顺 法 癌 涡 主 要 从 中 游 贺 柱 发 展 至 流 场 中 ,并 在 流 
场 中 发 生变 形 与 拉 伸 转变 为 其 他 方 回 的 旋涡 。 
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